Two-Leg Ladders and Carbon Nanotubes: Exact Properties at Finite Doping 
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Recently Lin, Balents, and Fisher have demonstrated that two-leg Hubbard ladders and armchair 
carbon nanotubes renormalize onto the integrable S0(8) Gross-Neveu model. We exploit this inte- 
grability to examine these systems in their doped phase. Using thermodynamic Bethe ansatz, we 
compute exactly both the spin and single particle gaps and the Luttinger parameter describing low 
energy excitations. We show both the spin gap and particle gap do not vanish at finite doping, 
while the Luttinger parameter remains close to its free fermionic value of 1. A similar set of con- 
clusions is drawn for the undoped systems' behaviour in a finite magnetic field. We also comment 
on the existence in these systems of the 7r-resonance, a hallmark of Zhang's SO (5) theory of high Tc 
superconductivity. 
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Carbon nanotubes are among a new class of quasi 
one-dimensional materials exciting tremendous interest 
in both the material science and physics communities |^] 
1^ [||. Nanotubes are built by wrapping graphite sheets 
into tubes of nanometer scale. They take a multitude 
of forms, ranging from crystalline ropes of nanotubes ar- 
ranged in triangular lattices to single multi-walled nan- 
otubes consisting of concentric cylinders to the simplest 
construction, a single- walled, isolated nanotube. The lat- 
ter promises to provide cleaner and longer realizations of 
quantum wires and as such promises to be the most rel- 
evant in new technologies 

Two-leg Hubbard ladders form another prominent 
class of quasi one-dimensional systems. Interest in these 
systems arises partly from the novel spin physics one 
expects in coupled chains and partly in that they 
provide an example of a Mott-insulator, a class of com- 
pounds that includes undoped high Tc cuprates. In con- 
trast to the latter, here long-range antiferromagnetic or- 
der cannot occur and there is a finite gap to spin ex- 
citations. As such they represent a spin liquid. Their 
basic doped behaviour has been examined both through 
weak coupling analytical treatments ^ and through 
Monte Carlo and density matrix renormalization group 
(DMRG) techniques d): the spin gap decreases in mag- 
nitude upon doping and the system exhibits quasi long 
range d-wave superconducting pairing correlations. We 
analyze this behaviour exactly in the following. 

Remarkably a certain class of the single walled nan- 
otubes behave as two- leg Hubbard ladders [||. Single 
walled tubes are classified by a superlattice translation 
vector, denoted by (N,M), specifying how the graphite 
sheet is to be rolled up in forming the tube. "Arm- 
chair" tubes, described by a superlattice vector (N,N), 
have been shown to be identical to a two-leg ladder with 
an effective Hubbard interaction of U/N. Our findings 
for doped ladders thus apply equally to doped armchair 
nanotubes. Due to their relative cleanness, it may well 



be in such systems that our results are most readily ob- 
served. Electronic properties of the tubes can be eas- 
ily measured by attaching metallic leads or by tunneling 
through an STM tip. The feasibility of the latter has 
been recently demonstrated by J. Wildoer et al. Ex- 
act results for the tunneling I{V)-cmve into an infinitely 
long tube have recently been obtained by two of us . 

Even more remarkable is the finding by Lin, Balents, 
and Fisher that the two- leg Hubbard ladder, and thus 
the armchair nanotubes, renormalize in weak coupling 
onto an SO (8) Gross-Neveu model. This equivalence is 
established by considering generic interactions. Under 
a one-loop renormalization group flow, the interacting 
Hamiltonian flows onto an SO (8) Gross-Neveu model, 

H= f (VI/1^Z9,*,L - ^in^^.:.^aR (1) 

^ a=l 



a theory of four interacting Dirac fermions, {ipRaiipLa)- 

SO (8) is an immense symmetry, usually found only 
in the domain of grand unified theories and supersym- 
metric strings In condensed matter systems we are 
usually limited to considering systems with charge con- 
servation (U(l) symmetry) and perhaps rotational invari- 
ance (SU(2) symmetry). But we are doubly blessed. The 
SO (8) symmetry resides in an integrable model. Often 
integrable models are considered to be of little use in de- 
scribing experiments in that they are highly tuned and 
thus are not expected to describe generic behavior in ex- 
perimental systems. However here we avoid this diffi- 
culty: the system is attracted to the integrable SO (8) 
Gross-Neveu model under a renormalization group flow. 
A similar situation was found in the tunneling between 
1^ = 1/3 quantum Hall edge states 

Intcgrability provides us with a tremendous amount 
of information. The key to organizing this information 
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is the exact knowledge of the basis of eigenstates of the 
fully interacting system. There is a precise notion of a 
"particle" in an integrable system akin to fermi liquid 
theory. However unlike fermi liquid theory the particle 
is perfectly stable. There is no particle decay or produc- 
tion; all scattering processes are encoded in two particle 
scattering matrices, S. 

Particles in the SO (8) Gross-Neveu model are orga- 
nized into three octets and a set of 29 bosonic particles. 
One octet is of fundamental (Majorana) fermions and 
is directly related to the particle/hole excitations associ- 
ated with the four Dirac fermions in eqn.(p. The remain- 
ing two octets are of fermionic "kink" particles. The 29 
bosons consist of a rank two anti-symmetric tensor with 
a single scalar boson, together forming a single represen- 
tation of the SO (8) Yangian. The masses of these parti- 
cle are generated dynamically through the Gross-Neveu 
interactions (the quartic term in cqn.(^). In terms of 
the ladder parameters, the mass scale, to, is governed 
by, m ~ te~*/'^, where t is the ladder bandwidth and U 
is the typical coupling strength. The fermions all have 
mass, TO, whereas the bosonic particles have mass, ^/3m. 

Knowledge of these masses together with knowledge 
of their quantum numbers was used to make a variety 
of predictions in with respect to the correlation func- 
tions at low temperature and zero doping. Largely qual- 
itative, these results were extended in ]lO| , where matrix 
elements were computed exactly in order to obtain exact 
forms for the correlators. Here we consider the doped 
system at zero temperature. Though more difficult theo- 
retically, it is more interesting experimentally. Beyond a 
critical value of the chemical potential, the gap to charge- 
two excitations goes to zero and the system becomes con- 
ducting. 

To model the doping, we add a chemical potential 
term, Hq = — /iQ, to the Gross-Neveu Hamiltonian 
where Q is the total electric charge. As Q is conserved 
by the original Hamiltonian, this addition does not spoil 
integrability. Rather it only lifts the degeneracy of par- 
ticles with differing charges. To determine the charge of 
the octet of fundamental (Majorana) fermions, we orga- 
nize them into four Dirac fermions (only in doing this do 
they have well defined charges). Of the four, one carries 
electric charge 2, the so-called cooperon, and one car- 
ries spin Sz = 1, the magnon. All four are associated 
with two-electron excitations. While the fundamental 
fermions carry but a single quantum number, the two 
octets of kink particles each carry ±1/2 of each of the 
four above quantum numbers. In particular, kinks carry 
charge, Q = ±1, and spin, Sz — ±1/2, and so a subset 
of them are identified with the physical electronic excita- 
tions on the ladder [Q. The 29 bosonic particles can be 
thought of as bound states of the fundamental fermions, 
and so carry two quantum numbers, for example, Q = 2 
and Sz = i- 

As the chemical potential is turned on, a particle with 



charge Q lowers its energy to to — /zQ. As the energy of 
a particular particle becomes negative, the ground state 
of the system changes. Once empty, it now begins to fill 
with particles, much as if the fermi energy was moved 
above the bottom of the conduction band in a semi- 
conductor. The first particle that will enter the ground 
state is the one with the highest charge to mass ratio. 
This is the cooperon. The ground state is pictured in 
Fig. 1 for 2/i > TO, where it consists solely of cooperons. 






bosonic bound states 
with Q=2 



FIG. 1. Schematic of the doped ground state. 

Seemingly as fi is increased further, the kinks and bosonic 
bound states will begin to fill into the ground state. How- 
ever the increase in chemical potential is counteracted 
by the repulsion felt by the kinks and bound states to 
the cooperons. It turns out that this repulsive energy is 
greater than the chemical potential for arbitrary fj,, and 
the ground state remains solely filled with cooperons, as 
follows from the general results of [Q. We emphasize 
that this is an exact result. 

We now proceed to characterize quantitatively the 
ground state through thermodynamic Bethe ansatz 
(TBA). At the root of the ansatz lies the insistence that 
the cooperon wave function be consistent with the scat- 
tering amplitudes. Sec (given in ||l^), of two cooperons: 

s r(i+f)r(i/2-f)r(5/6-f)r(i/3+f) 
> r(i-f )r(i/2+f )r(5/6+f )r(i/3-f )■ ^ ' 

In this way we can derive a density of occupied states 
for the cooperons, pc, and a corresponding energy, Cc, 
of the cooperons. We parameterize the non-interacting 
energy and momentum of a particle of mass, to, by 
E = TOC0sh(6'), p = TOsinh(0). With this parameteri- 
zation, the use of periodic boundary conditions results 
in an equation which determines the density of occupied 
states per unit length for \9\ < B: 



TO 

Pi&) = 7^ cosh(0) + / deh{o')Tcc{e - e') 



(3) 



Here Tec — -^dg In Sec is a measure of the strength of the 
cooperon-cooperon interaction; if Tec = we recover the 
density of occupied states for free fermions. ±B marks 
out the fermi energy. For \6\ > B, pc{0) — 0. One 
determines the interacting cooperon energy, ec(^), by as- 
sociating it with the occupied density of states at T 
via Pc = 1/(1 + exp{ec/T)). Demanding that the free en- 
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ergy be minimized and then taking T ^ 0, one obtains 
Ec in terms of the cooperon-cooperon S-matrix: 



de'Tr 



TOCOsh(e') - 2^. (4) 



-B 



Solving the constraint, £c{B) = 0, determines B. 

Having described the ground state, we now look at ex- 
citations over the ground state, that is the gaps. We 
have already seen that the gap to a charge-two excita- 
tion is zero; we can add a cooperon at an infinitesmial 
energy above the fermi energy. More interesting are the 
gaps to spin 1 and single particle excitations. As empha- 
sized above, these gaps never vanish at finite doping, a 
rigourous result. 

We first consider the spin gap. As. To compute this 
gap we want to identify the Q=0, Sz — I excitations. 
There are two possibilities: 1) adding a Q = 0, 5*2 = 1 
magnon; and 2) removing a cooperon from the sea and 
adding a bosonic bound-state with Q = 2 and = I. 
For 2fj, < m, the first excitation is the only one possible; 
the sea is empty. The spin gap is then — m. As we 
cross the threshold, 2/i — m, the Fermi sea begins to fill. 
At the threshold, the second process has a gap of only 
m{V^— 1), the energy of the cooperon hole added to the 
energy of the bound state 0. It is thus lower than the 
first process, and the spin gap is determined by this exci- 
tation. We can push beyond threshold and compute the 
gap for all values of the doping: previous studies 10 
have shown how to obtain the exact ground state of the 
theory at finite doping. The lowest Sz — I, Q — excita- 
tion is given by removing a cooperon at the fermi surface 
with rapidity. Oh = B, and adding a bound state with 
rapidity 9t, = 0. With these changes the particle density 
is modified to pc{0\e'' , 0'') . pc{e\e^\e^) is determined by 
an equation similar to (0), 



Pc(PWh, 



-B 



de'Tcc{e-e')pc{e'\eh,eb) = 
T,t,{e-et,)-T,,{0-0H). (5) 



Here we have used Trh — InS'^;,, where 



(6) 



Scb is the cooperon-bound state S-matrix jl^l Q. 

Having the new density, the energy of the excited state 
is then straightforward to compute: it is determined by 
the difference between the energy of the excited system 
and the energy of the original ground state. Given that 
the energy of the cooperon sea is L J d0pc{0) {cosh{0) — 
2p/m), the spin gap equals 



\/3 — cosh B 



m 



d0 5pc{0)[cosM 



2^), (7) 

TO 



where 5~pc{0) = Pci0\0h,Oi,) - Pc{0). 

The single particle gap, defined as the energy to add 
an electron with Q=l (i.e. a kink) to the system, follows 
from similar considerations. For 2p < to, the gap is that 
of adding a kink in the absence of a fermi sea, Af — m—p. 
Crossing the threshold, 2p > m, the sea begins to fill in 
and the interactions between the kink and the cooperons 
must be taken into account. The scattering between the 
cooperons and kinks, Sck is given by 



g,i|L,(g) + i^r(| + f)r(|-f) 
+ fL,(^)_^r(|-^)r(| + |^) 



Sck{Q) 



, (8) 



with s{0) = sinh(^?). The same procedure as before leads 
to the single particle gap, Af. 

-^ = 1-^+ d0 Spc{0){cosh{9)~2^). (9) 

TO TO J_g TO 

Here dpc is given by removing Tec and replacing Tcb by 
Fc/c = 9e In S'cfe/27ri, both on the r.h.s. of eqn. (5). 

Both of these gaps are shown in Fig. 2. The gaps have 
been plotted vs. D = / d0pc{0), the density of cooperons 
(that is, the carrier density). We observe that there is 
a uniform decrease with increasing doping, as expected 
for this system. For large D/m, the gaps approach zero 
asymptotically without crossing, as can be established 
analytically. 




FIG. 2. Plot of the spin and single particle gap. 

Although the gap to charge two excitations is zero, these 
excitations still encode interesting physics. The massless 
charge excitations occur near the fermi rapidity, — ±B. 
These excitations can be described by two chiral bosons, 
4'L/R.i (one for each fermi point). That the cooperon 
sea is interacting is reflected in the Luttinger parameter, 
K , that governs these massless bosonic excitations. The 
propagator of 4>l/r is given in terms of K. Focusing on 
(pL we have, {(j)L(z)4'L{w)) = —Khi[z—w). The operator 
creating a Cooper pair is thus Oc — e**^^/^, and so is 
governed by the propagator 



(10) 



where z = vpT -\- ix. K thus describes the long distance 
asymptotics of this correlator. 
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We can easily determine K in terms of the renormal- 
ized charge, Qr = — 9^ec|e=s/2 = Z, of the elementary 
excitations near the fermi surface. The basic idea, dis- 
cussed at length in is that the low energy excita- 
tions near the Fermi rapidity are just free fermions - that 
is, the dressed S matrix at low energy is just —1. This 
means that the only non vanishing matrix elements of 
the current operator, — dz4>L, in this basis involve 
particle- hole pairs, which, parameterizing the energy of 
the massless excitations by E{6) = e^, are given by. 



(11) 



The constant, c, is related to the interaction parame- 
ter, K, as we expect {jL{z)jL{0)) — ~K/z^. Using 
the matrix element to compute this correlator, we find 
K = —c^/Att'^. On the other hand, c is related via 
c = i2'!TQii to the renormalized charge, since integrating 
( |lT| ) over all x gives the value of Qr on the one particle 
states. One thus finds K = Qj^ = Z^. K, using eqn. 
(4), is then readily computed with the results plotted in 
Figure 3. With K < 1, the effective interactions are re- 
pulsive but because the cooperon-cooperon interaction is 
weak, the Luttinger parameter does not vary significantly 
from its free fermionic value of K = 1. Interestingly, the 
plot is marked by a dip at about D /ni ^ .4 before begin- 
ning an approach (albeit slowly) to an asymptotic value 
of = 1. 




FIG. 3. Plot of the Luttinger parameter. 

We point out that all of the above analysis applies 
equally to the system in a magnetic field, H, at zero dop- 
ing. Rather than coupling to the cooperon, the magnetic 
field couples to the magnon. However as the magnon 
and cooperon are related by an S0(8) rotation, all of the 
above equations hold under relabeling. In this scenario, 
the ground state consists of interacting magnons with 
zero spin gap. becomes the charge gap while re- 
mains the gap to single particle excitations (with = 1). 
Here K describes the low energy spin 1 excitations about 
the fermi surface. 

To conclude, we would like to make a few comments 
about the equivalent of the "tt resonance" discussed in 
Zhang p6l. Understanding singularities quantitatively 



is beyond the scope of this note; however, certain simple 
threshold effects can easily be analyzed. In particular, let 
us consider the spin spectral function. This can be com- 
puted, at least in principle, by inserting a complete set 
of states between the two spin operators, corresponding 
to various physical processes. We consider in particular 
spin excitations with energy Such excitations can be 
obtained by applying Zhang's 7r-operator to low energy 
Q — 2 excitations. The Q = 2 excitations are formed 
by adding a single cooperon to the ground state directly 
above the fermi sea together with cooperon particle-hole 
pairs. As there are an infinite number of such excitations, 
we expect the spin-spin correlator to exhibit, at momen- 
tum (7r,7r), a singularity of the type \uj — 2^P6(w — 2/z), 
akin to that found in (t). However unlike we expect 
this singularity to be robust: with integrability, mas- 
sive excitations well above the continuum can still be 
stable because they are protected by higher order con- 
served quantities preventing their decay, although clearly, 
non-integrable perturbations will give these excitations a 
small finite lifetime. 
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